Two component integrable systems modelling 
shallow water waves: the constant vorticity case 



Rossen 

School of Mathematical Sciences, 
Kevin Street, 



I. Ivano\0 

Dublin Institute of Technology, 
Dublin 8, Ireland 



Abstract 

In this contribution we describe the role of several two-component in- 
tegrable systems in the classical problem of shallow water waves. The 
starting point in our derivation is the Euler equation for an incompress- 
ible fluid, the equation of mass conservation, the simplest bottom and 
surface conditions and the constant vorticity condition. The approximate 
model equations are generated by introduction of suitable scalings and by 
truncating asymptotic expansions of the quantities to appropriate order. 
The so obtained equations can be related to three different integrable sys- 
tems: a two component generalization of the Camassa-Holm equation, the 
Zakharov-Ito system and the Kaup-Boussinesq system. 

The significance of the results is the inclusion of vorticity, an important 
feature of water waves that has been given increasing attention during 
the last decade. The presented investigation shows how - up to a certain 
order - the model equations relate to the shear flow upon which the wave 
resides. In particular, it shows exactly how the constant vorticity affects 
the equations. 
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1 Introduction 

The integrable nonlinear equations are used extensively as approximate models 
in hydrodynamics. They describe in a relatively simple way the competition be- 
tween nonlinear and dispersive effects. The best known example in this regard is 
the Korteweg-de Vries (KdV) equation. The use of term integrable corresponds 
to the idea that such equations are in some sense exactly solvable and exhibit 
global regular solutions. This feature is very important for applications where 
in general analytical results are preferable to numerical computations. 

The Camassa-Holm (CH) and Degasperis-Procesi (DP) equations [TJ^GI] are 
another two integrable equations with application in the theory of water waves 
[H HI El El E] • The excitement that greeted the CH and DP equations 
is due to their non-standard properties that set them apart from the classical 
soliton equations such as KdV. The first most remarkable of these properties is 
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the presence of multi-soliton solutions consisting of a train of peaked solitary 
waves (or 'peakons') ED [TO]. Another remarkable property of the CH and DP 
equations is the occurrence of breaking waves [TJ [TTJ CE3 EES EH] (i.e. a solution 
that remains bounded while its slope becomes unbounded in finite time [T5]) as 
well as that of smooth solutions defined for all times [TOl EEl [17] . 

In many recent publications the problem of water waves with nonzero vor- 
ticity and especially with constant vorticity is under investigation - e.g. see the 
publications EE3 EEi EM E3 EM EM EM Ell [Ml W\ and the references therein. 
The nonzero vorticity case arises for example in situations with underlying shear 
flow [6j. 

Our aim is to describe the derivation of shallow water model equations for the 
constant vorticity case and to demonstrate how these equations can be related to 
some other integrable systems: a two component generalization of the Camassa- 
Holm equation [28] , Zakharov-Ito system EM [3D] and Kaup-Boussinesq system 
[31 j . A starting point in our derivation are the equations that express the 
constant vorticity and the mass conservation. Another approach, based on the 
Green-Naghdi approximation is used in an alternative derivation of the two 
component Camassa-Holm equation in [32] . where the occurrence of solutions 
in the form of breaking waves is also established. The Kaup-Boussinesq system 
is used as a model in hydrodynamics under slightly different assumptions in 
[Ill[3ll[33]. 



2 Governing equations for the inviscid fluid mo- 
tion 

The motion of inviscid fluid with a constant density p is described by the Eulcr's 
equations: 

<9v 1 

w + (v.V)v = --VP + g, 

V-v = 0, 

where v(x,y, z,t) is the velocity of the fluid at the point (x,y,z) at the time 
t, P(x, y, z, t) is the pressure in the fluid, g = (0, 0, — g) is the constant Earth's 
gravity acceleration. 

Consider now a motion of a shallow water over a flat bottom, which is located 
at z = (Fig. 1). We assume that the motion is in the x-direction, and that 
the physical variables do not depend on y. 




Figure 1: Water waves: general notations 
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Let h be the mean level of the water and let r)(x,t) describes the shape of 
the water surface, i.e. the deviation from the average level. The pressure is 
P(x, z, t) = Pa + pg(h — z) + p(x, z, t), where Pa is the constant atmospheric 
pressure, and p is a pressure variable, measuring the deviation from the hydro- 
static pressure distribution. 

On the surface z = h+rj, P — Pa and therefore p = rjpg. Taking v = (u, 0, w) 
we can write the kinematic condition on the surface as w — r\t + ur\ x on z = h + r/ 
[34] . Finally, there is no vertical velocity at the bottom, thus w = on z = 0. 
All these equations can be written as a system 

1 

u t + uu x + wu z = —-p x , 

p 
1 

w t + uw x + ww z = — p z , 

p 

u x + w z = 0, 

w = r]t+ur) x , p = ripg on z = h + rj, 
w = on z = 0. 

Let us introduce now dimensionless parameters e — a/h and 8 = h/X, where 
a is the typical amplitude of the wave and A is the typical wavelength of the wave. 
Now we can introduce dimensionless quantities, according to the magnitude of 
the physical quantities, see [5J [Ml [3S] for details: x — > Xx, z — > zh, t — > -^^t, 

rj — > ar], u — > e^/ghu, w — > ed^/ghw, p — > spgh. This scaling is due to the 
observation that both w and p are proportional to e i.e. the wave amplitude, 
since at undisturbed surface (e = 0) both w — and p = 0. The system in the 
new, dimensionless variables is 



u t + e(uu x + wu z ) = -p x , 

6 2 ( w t + e(uw x + ww z )) = -p z , 

u x + w z = 0, 

w — r)t+sur)x, p — i] on z = l + er), 

w — on z = 0. 



3 Waves in the presence of shear 

So far no assumptions have been made on the presence of shear. Now let us 
notice that there is an exact solution of the governing equations of the form 
u = U (z), < z < h, w = 0, p = 0, rj = 0. This solution represents an arbitrary 
underlying 'shear' flow [6]. Let us consider waves in the presence of a shear flow. 
In such case the horizontal velocity of the flow will be U{z) + u. The scaling 

for such solution is clearly u — > y/gh(jj(z) + euj , where u on the left-hand side 
is the horizontal velocity before the initial scaling, and the scaling for the other 
variables is as before. Thus, in this case we have (the prime denotes derivative 
with respect to z): 

u t + Uu x + wU' + e(uu x + wu z ) = ~p Xl (1) 
S 2 ( w t + Uw x + e(uw x + ww z )) = -p z , (2) 
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u x + w z = 0, (3) 
w = f)t + [U + eu)rj Xl p = i) on z = l + et) (4) 
w = on 2 = 0. (5) 

We consider only the simplest nontrivial case: a linear shear, U(z) = Az, 
where A is a constant ( < z < 1). We choose A > 0, so that the underlying flow 
is propagating in the positive direction of the x-coordinate. Burns condition [36] 
gives the following expression for c, the speed of the travelling waves in linear 
approximation : 

c = \(A± ^4 + A 2 ). (6) 

This expression will be derived again in our further considerations, e.g. see (ITS]) . 
Note that if there is no shear (A = 0), then c = ±1. 

Before the scaling the vorticity is u> — (U + u) z — w x or in terms of the 
rescaled variables {uo — > \fh/g uo), 

uj = A + e(u z — S 2 uu x ). 

We are looking for a solution with constant vorticity uo = A, and therefore we 
require that 

Uz - S 2 w x = 0. (7) 

This assumption amounts to considering approximate wave-solutions that 
are interactions of an underlying shear flow and an irrotational disturbance 
thereof. From ([7]), and © we obtain 

u = u -5 2 ^u 0xx +O(e 2 ,5\eS 2 ) : (8) 
z 3 

w = -zu ax + S 2 —u 0xxx + O(e 2 ,5 4 ,eS 2 ), (9) 
6 

where uo(x,t) is the leading order approximation for u. Note that uq does not 
depend on z since from ([7]) it follows that u z — when S — > 0. 
From (JU) with © and © we obtain 



i] t + Ar] x 



A 
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ignoring terms of order C(e 2 , S 4 , eS 2 ). 
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(1 + er])u {) + e-rf -S 2 -u 0xxx = 0, (10) 



G 



From ([H, ([S]) and ^ we have (again, ignoring terms of order 0(e 2 ,5 4 ,ed 2 )) 



P = V 



ri - z z l - z A 

uoxt H z — ^ u 02 



2 3 
Then ([T]) gives (note that there is no z-dependence!) 

/ 1 \ A 

[u - S 2 -u Qxx j^ + eu u Qx + r) x - 5 2 —u 0xxx = 0. (11) 

Letting both the parameters e and S to 0, we obtain from (fTU|) . (fTTj) the system 
of linear equations 

u ot + Vx = 0, (12) 
rjt + Ai] x + u 0x = 0, (13) 
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giving 



rj tt + Ar] tx - Vxx = 0. 



(14) 



The linear equation (|14p has a travelling wave solution r/ = r)(x — ct) with a 
velocity c satislying 



Ac- 1 = 0. 



(15) 



This gives the same solution for c that follows from the Burns condition ©. 
There is one positive and one negative solution, representing left and right 
running waves. We assume that we have only one of these waves, then (e.g. 
from P}) 

r 1 = cu Q + 0(s,6 2 ). (16) 



Let us introduce a new variable 

p = 1 + ear) + s 2 (3rf + sS 2 ^u 0xx , 



(17) 



for some constants a, 8 and 7. These constants will be determined in our 
further considerations. The variable p will be used instead of 77 as a tool for 
mathematical simplification of our equations. The expansion of p 2 in the same 
order of e and 5 2 is 



1 + e(2a)r] + s 2 (a 2 + 2/3)r; 2 + e5 2 {2 1 )u 0x 



(18) 



With this definition one can express rj in terms of p and write equation (110[) 
in the form (keeping only terms of order 0(e,S 2 )): 



Pt + Ap x 



as 



+ 5 2 {-{c- A)- -)u 0xa 
A 



(3 



(1 + eri)uo + e— r] + e—cu 
2 a 



= 0. 



One can eliminate the Moa^-term by choosing 
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Equation (|T9|) becomes 
p t + Ap x 



a 6(c-A)' 



Ac 3. , 

1 + e(l + + -)»? «o 
2 a 



0. 



Ac 3 
-2 + a 



With the choice 

a = 1 

we can write (|21j) in the form 

+ Ap x + ae(puo)a; = 0, 
which contains only the variables p and uq but not 77. 



(19) 



(20) 



(21) 



(22) 



(23) 
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4 Two component Camassa-Holm system 



In this section we proceed with our derivation in a direction that leads to a two 
component Camassa-Holm system. Expressing r\ in terms of p in (jlip we obtain 
(matching terms of order 0(e, S 2 )): 



m t + Am x - Au 0x + S 2 ( — + k(A - c) - -)u 0xxx 

V a/ 

+ e {i^^±li c A UoUOx + J-^) x = 0i (24) 
V a / lea 



where m — uq ~ 5 2 {\ + k)uq xx , k is arbitrary: we are adding and subtracting 
5 2 KUQ XX t, making use of (jTHJ) - Fixing 



1 / 7 A 



A — c \a 6 

leads to the disappearance of the uq xxx - term . 
The relations {20]) and {25]) give 



(25) 



K= 6^A){ A -7h)' (26) 

Thus equation (|24p can be written as (matching only terms up to the order 
0(e,S 2 )) 

m t + Am x - Au 0x + e\ ( 1 - - +2 ^ c 2 ) [2mu 0x + u m x ] + = 0. (27) 

Recall that m = Uo — 5 2 Buq XXi where (see (|2l)]) and {H}) 

1 _ A 2 -c 2 + 2 _ 1 

B ~ K+ 2~ 3(c-A) 2 ~3c 2 (c-A) 2 - (28) 



Note that B is always positive and the denominator in ([2"5]) nonzero (since c ^ A 
- see (JTSJ). The rescaling u -> sjUo, £ -> -75 i,t-» -7=* in {27]) and {23]) is 
now only for the sake of mathematical clarity and simplicity and gives: 



1/ a 2 + 2/3 ,\ . . 

m t + Am x - Au 0x + — 1 c 2mucta + t^m,! + pp x = 0, 

da V a / 

m = U - tiQxx; 

p t + Apa; + (pUo)^ = 0. 

Finally, we choose 
and thus 

m t + Am^ - Auo-c + 2mu 0x + u m x + pp x = 0, m = u —u xx, (30) 
Pt + Ap x + (pu ) x = 0. (31) 
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The constants a, (3 and 7 can be determined from the constraints (f2"2"l) . f2H) 
and Jig): 
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1 

3(1 + c 2 ) + 3 ' 

1 3 + c 2 



3(1 + c 2 ) 6 
1 



a, 



6(c-A) 



(32) 
(33) 
(34) 



Note that from (1321) it follows that a is always positive. 

Now let us express the original variable r\ in terms of the 'auxiliary' variable 
p. Before the rescaling we had aerj = p — 1 — e 2 f3c 2 u 2 ] — eo~ 2 ^uq xx . Since in 



the leading order 77 
rescaled variables 



cuq the rescaling of rj is rj 



thus in terms of the 



V = P 



1 



2 Tjl 

—^-u Q - B—Uo xx . 



With a Galilean transformation (that we use only to simplify our equations 
and to bring them to the form that is widely used), such that dt> = dt + Ad x , 
d X ' — 9 X (%' = x — At, t' — t) we obtain 

m t , - Au Qx > + 2mu x' + u m x > + pp x > = 0, m = u Q - u 0x 'x> (35) 
p t > + {pu ) x > = 0. (36) 

The system (|36p is an integrable 2-component Camassa-Holm system 

that appears in [28j , generalizing the famous Camassa-Holm equation PQ . 
Let us drop the primes and write it in the form 

m t - Au 0x + 2mu Qx + u m x + pp x = 0, m = u - u 0x x (37) 
Pt + {pu o ) x = 0. (38) 

It generalizes the Camassa-Holm equation [1 in a sense that it can be obtained 
from it via the obvious reduction p = 0. The system is integrable, since it can 
be written as a compatibility condition of two linear systems (Lax pair) with a 
spectral parameter £: 



*x*= ( -( 2 p 2 + ({m- 4) + 



1 



*i = 



1 

2C 



1 



The system is also bi-Hamiltonian. The first Poisson bracket is 

f [5Fi 5F? SFi SFo SFi SF 9 

{F lt F 2 } = - / i±{-Ad + md + dmf-p + -^pd-^ + -j^dp -^ 
J I dm dm dm dp dp dm 

for the Hamiltonian Hi = i f (uo{m — 4) + p 2 )dx. 



dx 



The second Poisson bracket is 
{Fi, F 2 } 2 = ■ 



— (d-d iS )— + —d— 
dm dm dp dp 



dx 
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for the Hamiltonian H2 = \ j {uop 2 + Uq + uqu^ x — Aufydx. It has two Casimirs: 
J pdx and J mdx. 

The system has an interesting interpretation in group-theoretical context. 
The first Poisson bracket gives rise to a Lie-algebraic structure. This fact is well 
studied in the case p = when the system coinsides with the Camassa-Holm 
equation 38, 39, 40JHT122]- Then the corresponding Lie algebra is the Virasoro 
algebra. 

By considering the expansions 

we obtain the following Lie-algebra with respect of the first Poisson bracket: 



i{L n ,L k } = (n - k)L n+k - 2nAn5 n+k , (39) 
i{p ni L k } = np n+k , (40) 
i{p n ,Ph} = (41) 

The Lie algebra (f3"9"]) - (|4"Tj) is a semidirect product of the Virasoro alge- 
bra (trit) (|39p with a central charge proportional to A, and the abelian algebra 
C°° (R) (jHJ) [37] . Note that the central extension of the Virasoro algebra con- 
tains only An term but not An 3 term since the Hamilton operator contains only 
the first derivative Ad. The system (|37[) . (|38[) represents the equations of the 
geodesic motion on the corresponding Lie group QJir x C°°(R) for the metric, 
defined by Hi: 



ll(«o.p)lr = / « + u L + P )dx, 

which is right-invariant under the natural group action. 



5 Zakharov-Ito system 

In this section we describe a derivation that matches the approximate equations 
(fTTjl . ([23| (where p is given in (fT7| ) to the integrable Zahkarov-Ito system [29| 
1301 SUS]: 



Mot - 4:ku 0x + uo X xx + 3u uox + PPx = 0. (42) 
Pi + {u Q p)x = 0, (43) 

where k is an arbitrary constant. The system is formally integrable by the virtue 
of the Lax pair [13 02] 

= (c- y + fc-r 1 ^)*, *t = -(4c+«o)*x + 5«ito*. 

Equation (fTT|) can be written in the form (cf. (fl6|) and (fT2|) ) 



8 



u ot + S 2 Ku 0xxx + e(l - +2 ^uqUox + — = 0. (44) 
V a / as 

where if is a constant, given by [see (|20|) and (fT5|) ] 

^=f-4~ 2 -^ c - A )- ( 45 ) 
2 3 a 3 

For one of the roots of the equation in ((6|) , K is positive and for the other it 
is negative. We can fix the constants a, (3 and 7 by the conditions (f2"2")) and 

i(i-^)=., 

da V a / 

which are formally the same as those for the Camassa-Holm system, giving the 
same expressions (j3"2")l - (j3~4"|) . 

The rescaling u — > tjjMo, £ — ► * — > m 631) an d (HSJ) gives: 

+ 3uou 0a ; + ppx = 0, 

Pi + Ap x + (pu ) x = 0. 

Note that a coordinate change (x,t) — > i(x,t) maps into another system with 
real variables and therefore if K < we still can apply formally the above 
rescaling. One can use a Galilean transformation x' = x — At, t' = t to obtain 



uot' - Au 0x ' + u 0x , x > x , + 3u u 0x ' + pp x > = 0, 
pt> + (pu a ) x > = 0, 

that matches the Zakharov-Ito system B2I , (|4"3")) if the constant is chosen to 
be k = A/4. The 'physical' variable 77 in terms of the 'auxiliary' p and uq (for 
the rescaled variables) is 

1 ^ 2 ^7 
V = p - 1 5-^0 - ^Oaa- 

6 Kaup-Boussinesq system 

Another integrable system matching the water waves asymptotic equations to 
the first order of the small parameters e, 8 is the Kaup-Boussinesq system [3TJ 
133] . In this section we describe briefly its derivation. Introducing 



2/1 A\ . 2 /l 



the equation (jTTj) can be written as 

Vt + sVV x + t] x = 
Equation (flU)) in the first order in e, <5 2 is 



1 



(46) 



Vt 



A 

Arj + (1 + erj)u + £-^V 2 



,1 



<5 -^uq xxx — 
6 
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and with a shift 77 — s- 77 — -i- it becomes 

T} t + e(l + 4~)( ? ? u o)x - S 2 ^u 0xxx = 0, 
2 

or 

r/t + e^^n.-^^ = 0. (47) 
2 

Further rescaling in ([4!^ and (j4"7f leads to the Kaup-Boussinesq system 

+ vv x + vx = o, m- \v xxx + l^L(riV) x = 0, 

which is integrable iff A — (c 2 = 1) with a Lax pair 

*t = -(C + ^0** + ^**- 

The integrability of the system, as well as the Inverse Scattering Method for 
it has been investigated firstly by D.J. Kaup (21] ■ His motivation has been to 
derive an integrable water-wave system with a second-order eigenvalue problem, 
which is readily solvable in comparison to the third-order eigenvalue problem 
for the Boussinesq equation. In our context, however, this system is relevant 
only in the case with zero vorticity. 

7 Discussion 

Apparently the described method can be used for other two-component inte- 
grable systems with a similar structure, e.g. see the classification in [46j . It is 
interesting to investigate further which specific properties of the original gov- 
erning equations are preserved in the 'integrable' approximate models. For 
example the 2-component Camassa-Holm system for certain initial data admits 
wave breaking [32 . Peakons do not occur in the case A — J3_2j and most 
certainly not in the case A ^ 0, due to the term with linear dispersion Auq x - 
However, in the 'short-wave limit' where m = —uq xx and A = peakon solu- 
tions are possible [32] . Recently a similar system with peakon solutions have 
been constructed in [IT] . 

The case with — pp x term (instead of +pp x ) in (J37J is also integrable [48 
and it is studied in [ 191 . The two component Camassa-Holm system appears 
also in plasma theory models [50[ 15 lj and in the theory of metamorphosis [52j . 
Other integrable multi-component generalizations of the Camassa-Holm equa- 
tion (including other two-component ones) are constructed in [48] ■ 
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